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1. INTRODUCTION 


It has been hypothesized that boundary pixels, the so-called iiiixels, are a 
major source of classification error in the various clustering and classifi- 
cation algorithms applied to LANDSAT data. This classification error is due 
to (1) the distortion of the statistics for the classes identified by the 
algorithm caused by the inclusion of different targets and (H) because the 
label assignment for nonhomogeneous areas is not well defined. It is expected 
that if the boundary pixels were removed, the accuracy of clustering and 
classification would be greatly improved. 

This document describes programs which generate an image file that has all 
ground truth boundary pixel spectral values set to one value. This image, 
when processed by LARSYS routines, gives classification and clustering maps 
with all boundary pixels assigned to one class. 

Using these programs the performance of clustering and classification 
procedures for pure pixels can bo tested. 

2. DI-SCIUPTION OF THE SOFTWARE 

The ground truth data are available at JSC, thus the ground truth processing 
takes place on the image processor in the Data Techniques Laboratory. The 
spectral data are available both at JSC and LARS. Since the clustering and 
classification system is included in LARSYS at Purdue, the ground truth purity 
data have been merged with the spectral data at LARS. The flow of data is 
shown in figure 1. 

2.1 Image Processor Programs 

2.1.1 GROUND TRUTH INPUT 

The ground truth input comes from disk files installed by accuracy assessment 
software on a disk mounted on the second disk drive (DB2) of the image 
processor. This input is a digital map generated from ground truth data with 
six pixels for each LANDSAT pixel. These data are documented in "Format 
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Specifications for LACIE (Phase III) and Accuracy Accessinent Computer Data 
Products," LEMSCO- 12507. 

2.1.2 PIXEL PURITY IMAGE PROGRAM (PUROUT) 

2. 1.2.1 Linka ge 

PUROUT calls subroutine PURE at the entry points PURE, PUREl, PURE2 and ROLL. 
PUROUT calls subroutine STRIP at the entry points STIN and STRIP. 

2. 1.2. 2 Interface 

PUROUT communicates with PURE and STRIP through the common block /PURE/. 

2. 1.2. 3 Input 

PUROUT reads FILNAM.PAT to get the name of the imago data file on logical unit 2. 
PUROUT reads the image file for example DB2:[111,3]013579999.GT0, on Logical 
unit 1. 


2. 1.2. 4 Output 

PUROUT writes the image file of type PTl on logical unit 3. For example: 
C111,3]013579999.PT1). 

2. 1.2. 5 Storage 

Total space allocated 3184. 

2. 1.2. 6 Description 

PUROUT first reads the input file name and checks it interactively with the 
operator. Then the header is copied, unchanged, to the output file and 
processing begins. 

PUROUT holds three lines of data in the array BUF (392,9). Since each pixel 
corresponds to 6 subpixels each line occupies a 392 x 3 block of space. The 
data are read into sublines 7, 8, and 9. Subroutine PURE determines whether 
the pixels on the input line are pure. Pure means that all subpixels are in 
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the same class. Subroutlr.'e ROLL moves all the data up one line. Subroutine 
PUREl checks pure pixels in 4, 5, and 6 and marks those surrounded by 
subpixels of the same class as "more-pure". PURE2 checks pure pixels in 4, 5, 
and 6 and marks those surrounded by pure pixels of the same class as "most- 
pure". STRIP removes strip-fallow classes. As data are rolled to the top of 
the array it is written to the output file. 

2.1.3 SUBROUTINL PURE 

2 . 1 . 3 . 1 Linkage 

PURE has four entry points, PURE, PUREl, PURE2, and ROLL. 

2. 1.3. 2 Interface 

All information is transferred through the common block/PURE/. This block 
contains the input Byte data, array, BUF (392,9), and the output integer *2 
array LAB (196,3). 

2. 1.3.3 Input 
No input 


2. 1.3. 4 Output 
No output 


2. 1.3. 5 Stor»ge 

Total space allocated 2729 

2. 1.3.6 Description 

1. Entry PURE ~ A pixel is pure only if all subpixels are of the same 

subclass. Subroutine PURE marks pure pixels with "1" and impure pixels 
with "0". The majority label is also saved. PURE works on the bottom 
line, sublines 7, 8, and 9 of the input array. 
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2. Entry PUREl ~ PlIREl checks pure pixels in line 2 and sublines 4, 5, and 
6, The purity label is changed to "2" if all neighboring subpixels have 
the same label as the pixel. The figure below shows the order of 
checking. 
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3. Entry PURL2 ~ PURER cliecks more pure (Lab, el - "2") pixels in line 2 and 
sublines 4, b, and 6, and changes the label to "3" if all adjacent pixels 
have the same label. Pixels are checked In the order shown in the figure 
bel ow. 


1 

5 

3 

7 


8 

4 

6 

2 


4. Entry ROLL ROLL moves the dfrta up one line and three sublines, in 

preparation for new input. Line 1 and sublines 1, 2, and 3 are destroyed 
in the process. Sublines 7, B, and 9 are not cleared. 

2.1.4 SUBROUTINE STRIP 

2. 1.4.1 Linkage 

STRIP has two entry points, STRIP and STIN. 

2. 1.4. 2 Interface 


STRIP transfers data through the common block /PURE/ and through the calling 
arguments. 





2. 1.4. 3 Input 

Entry STIN reads the array ZAP (256) on logical unit 8 from the file specified 
in the array CRDFIL (32). 

2. 1.4. 4 Output 

STIN may type an error message. 

2. 1.4. 5 Storage 

Total space allocated 2778 

2. 1.4.6 Description 

1. STRIP (KK) 

STRIP changes the purity class of marked classes to "0". The marked 
classes are typically problem fields, strip fields, and non-inventoried 
fields. These fields can be marked in two different ways. Three 
different sets of class identifiers may be coded into the array TEST (8, 3), 
or the array ZAP (256) may be read in. If KK = 10, the ZAP alternative is 
used. If KK = 1, 2, or 3 the array TEST (*, KK) is used. 

2. STIN (CRDFIL) 

STIN reads the array ZAP from the file specified by CRDFIL. if ZAP (N) - 
0, class N is accepted, if ZAP (N) = 1, class N is marked and purity will 
be set to "0". 

2.1.5 TAPE GENERATION (TAPEOUT) 

The program TAPEOUT outputs a universal format tape from a disk file 
(reference Action Document, 63-3107-4402-16). 

2.1.6 TASKBUILDER COMMAND FILE 

The file PUROUT.CMD contains the taskbuilder commands needed to construct the 
PUROUT.TSK file. 


6 



Z.Z LARS PR O GMM S 

2.2.1 TAPC TRANSFER (TAPTRAN) 

TAPTRAN is a proijrain written and maintained by Purdue. It is documented in 
"LARS DATA - iOU Operator's Manual." 

2.2.2 PURITY IMAGE TAPE TO DISK (TPURlO) 

2. 2. 2.1 Linkage 
None 

2. 2. 2. 2 Interface 
None 

2. 2. 2. 3 Input 

TPURCO reads a universal format* 1 channel tape from unit 11. The line size 
is 90 rNTEGER*4 words or 360 Bytes. 

2. 2. 2. 4 Output 

The header is copied to unit 13. The data are copied to unit 12. A small 
report is written on unit 6. 

2 . 2 . 2 . 5 Storage 
Program Size =5 4050. 

2. 2. 2. 6 Description 

TPURCO reads the tape and copies it to disk files. 

2.2.3 SPECTRAL VALUE TAPE TO DISK (TAPCOP) 

2. 2. 3.1 Linkage 

TAPCOP calls GETACQ, RTEERR, and TOPRD. These are all Purdue maintained 
routines. Documentation can be found in LARS Program abstract 11 for module 
TAPOP and LARS Program Abstract 2020 for module GTINFO. 
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Z,Z>Z,t Interface 


Interface is through the calling arguments and the tape mounted by GETACQ. 
^•2*3.3 Input 

TAPCOP interactively gets the segment name and date and TOPRD reads the tape 
mounted by GETACQ, 

2. 2. 3. 4 Output 

The header of the universal format input tape is written on unit 13. The 4- 
channel spectral data are written on unit 12. 

2.2.3. 5 Storage 
Program size - 4822. 

2. 2. 3. 6 Description 

TAPCOP interactively gets segment and data. These are passed to GETACQ which 
mounts the correct LARS library tape and positions it at the correct file. 
RTEERR decodes the error flag returned by GETACQ. If there is no error, TOPRD 
reads the tape which is then written to a disk file for further processing. 

2.2.4 FILE MERGE (DSKRED) 

2. 2. 4.1 Linkage 
DSKRED calls TRNSLT. 

2. 2. 4. 2 Interface 

The interface is through the calling arguments. 

2. 2. 4. 3 Input 

DSKRED interactively gets the desired purity class from the terminal. A 
universal header is read from unit 17. Four-channel spectral data are read 
from unit 18 and 1-channel purity data are read from unit 19. 
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2.2»4.4 Output 

ea— 

A universal format image tape file is written on unit 20. 


2. 2. 4. 5 Storage 
Program size ^ 9910, 

2. 2. 4. 6 Description 

DSKRfl) reads the spectral values and purity values. For those pixels with 
less purity than is desired, the spectral values are changed to ch 1 « 0, ch 2 
» 0, ch 3 0, and ch 4 255, 

Then these data are written out in universal format. 

2.2.5 BYTE MANIPULATION (TRNSLT) 

2. 2. 5.1 Linkage 

Subroutine TRNSLT does not call any other program. 

2. 2. 5. 2 Interface 

All data are passed through the calling arguments. 

TRNSLT (DUF, PUF, DHl, DH2, DH3, D114, PH, OPTION) 

DUF (225) spectral value input line 
PUF (90) purity value input line 
DHl (196) - DH4 (196) output values 
PH (196) purity output values 
OPTION. If OPTION ^ 1 DH 

is KAUTH transformed 

If OPTION - 2 DH 

is LANDSAT 3 corrected 

and KAUTH transformed 

If OPTION - 0 DH is raw channel values. 
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None 


2. 2. 5. 4 Output 
None 

2. 2. 5. 5 Storage 
Program size - 0202. 

2. 2. 5. 6 Description 

IIUWX'-TTK*!— *ie:TCKSKsm-jK!mt» 

Subroutine TRNSLT converts one line of spectral data in bytes to four-integer 
arrays, and also one line of purity data in bytes to an integer array. The 
spectral output may be raw channel values, KAUTIl transformed values or 
Lockheed/PMSCO's LANDSAF 3 corrected KAUTIl transformed values. 

The data are placed in LOGICAL*! arrays by equivalence statements and then 
assnsin?<:i to integer arrays. 

2.:-. 6 EXECUTIVL ROUTINES 

For the programs TPURCO TAPCOP and DSKRED there are EXEC files with the same 
names which give the required FILEDEF commands and start execution. In 
addition RTF EXEC may be executed to give the needed GETDISK commands. 
Subroutines GLTACQ, RTEERR, and TOPRD reside on JSC19A . 

2-3 MODIFICATIONS TO LARSYS ROUTINES 

To run LARSYS on the output file of DSKRED the supervisor program, MONTOR, is 
used. Also, three of the LARSYS subroutines required slight modification. 
These modified programs reside on 0SC8O8. 
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2.3.1 WRTHED 


The information saved by LARSYS for the output tape header was deemed 
inadequate. Therefore code was added to read the input file header and write 
it to the output file after two small changes were made. 

WRTHED calls the new subroutine IDTE for current date. 

2.3.2 PSPPAT 

Because the channel 4 value for impure pixels is set to 255, the accumulator 
register for sum of channel 4 squared has excessive error. It was necessary 
to change this variable to REAL*8, double precision to avoid excessive error. 

2.3.3 COVPAT 

Because the values for all impure pixels are the same, the covariance matrix 
for impure pixels is singular. LARSYS rejects singular covariance matrices. 
Code was modified in subroutine COVPAT to insert a nonsingular covariance 
matrix whenever mean channel 4 exceeds 250 counts. 

2.3.4 CURRENT DATE (IDTE) 

2. 3. 4.1, Li nkage 

IDTE calls GTDATL' and I MONTH. 

2. 3. 4. 2 Interface 

All interface is through the calling arguments. 

NA is month number 
NR is day of month 
NC is year. 

For the system routine GTDATE the array data are in A format. For example, 
printing DATE as 3A4 gives June 14, 1980. 

2. 3. 4. 3 Input 
None. 
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g.3»4.4 

None. 


3. 3. 4. 5 Sto rnae 
PrmjraiH Size SP4. 

P.3.4.6 Description 

■: *.r* rrv tastSi 

Subroutine IDll obtains the currrni date in Alphanumeric format. Conversion 
of the numbers is done by writing and rereading]. IMONIil is called to convert 
the month. 

2.3.B MONTH CONVI.RSION (IMONTll) 

2. 3.5.1 Linkage 
None, 

?.3v5.2 I nterf ace 

Date is a 4-cl<nracter month name. I is the inteyer month number. 

2. 3. 5.3 Inmit 

ii rr iwt ' . 

None . 

2.3. 5.4 Output 

Possible error statement ott unit 6. 

2. 3. 5.5 Storgye 
Program size - 495. 

2.3. 5.6 Oescriplion 

IMONTll compares the month name i.o tost values until n match is found. If no 
match is found an error statement is generated. 
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